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ABSTRACT 

We here investigate the possibility that the observed collimated outflows in 
luminous infrared galaxies (LIGs) and some Seyfert galaxies can be produced in 
a starburst (SB) environment. In the former source class, in particular, there 
seems to be some observational evidence for the presence of nuclear SBs in some 
objects (e.g. Scoville, Lonsdale, & Lonsdale 1998). 

A nuclear disk can be quickly produced by gas infall during star formation 
in a rotating, stellar cluster. We find that massive nuclear SBs with core 
disk masses M d ~ 10 8 — 10 9 M Q , and supernova rates u S n — 5 x 1CT 3 — 2 
yr _1 (which are consistent with the vsn values inferred from the observed 
non-thermal radio power in source candidates) may inject kinetic energies which 
are high enough to blow out directed flows from the accreting disk surface, 
within the SB lifetimes. In our models, the acceleration and collimation of 
the nuclear outflow are provided by magnetic fields anchored into the rotating 
SB-disk. The emerging outflow carries a kinetic power that is only a small 
fraction (a few percent) of the supernovae energy rate produced in the SB. 
Based on conditions determined from observed outflows and disks, we find that 
moderate disk magnetic fields (^8x 1CT 4 G) are able to accelerate the outflows 
up to the observed terminal velocities (^ few 100 km s -1 in the case of the 
Seyfert galaxies, and ~ 400 - 950 km s -1 in the case of the LIGs). The outflow 
is produced within a wind zone in the disk of radius <^ 100 pc in the the LIGs, 
and <^ 10 pc in the Seyferts, with wind mass loss to disk accretion rate ratios 
M w /Md 'k, 0.1 (where M d ~ 100 M yr _1 ). The observation of rotating nuclear 
disks of gas within ~ few 100 pc scales, like that in the prototype LIG, Arp 220 
(for which M d ~ 100 M yr~ x ), and magnetized outflows in Sey galaxies and 
LIGs (with terminal B ~ 10~ 5 G at the kpc-scales), provide some observational 
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support for the magneto-centrifugal disk picture drawn here. 

Subject headings: galaxies: active - ISM: jets and outflows - galaxies: 
starburst, Seyfert galaxies, LIGs 
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1. Introduction 

Active Galactic Nuclei (AGN) emit large amounts of energy (10 43-48 ergs s -1 , from 
Seyfert galaxies to QSOs) over a large range of the spectrum coming from a small volume. 
The nature of the central engine that produces these extraordinary amounts of energy 
remains virtually unknown. Current AGN theories include: i) black holes (BHs) accreting 
matter from the surroundings (e.g. Rees 1984); ii) Starburst (SB) models in which activity 
is due to star formation in a central cluster (e.g. Terlevich et al. 1992 - who have revived old 
theories, e.g., of Shklovskii 1960, Field 1964); and iii) BHs surrounded by dense starbursts 
(e.g. Perry & Dyson 1985, Scoville & Norman 1988, Williams & Perry 1994). 

While BH models (with masses ~ 10 6 — 1O 1O M ) accreting matter from the 
surroundings have the power to produce the observed luminosities and naturally provide 
a preferred axis for collimating jets, they seem to explain unsatisfactorily, e.g., the broad 
emission-line regions (BELRs) of AGN (as they require a fine-tuning to fit the observations). 
In Starburst models, on the other hand, strongly radiative compact supernova remnants 
(cSNRs) are the source of the broad emission lines. This model has been proved plausible 
in explaining the observed emission line spectrum, the UV-optical continuum and related 
variability (e.g., Terlevich et al. 1992), but is difficult to be applied to the whole population 
of AGN. In particular, they do not provide a natural mechanism for jet production and 
collimation. We here address this question. 

It seems plausible that active star formation and a compact AGN activity may coexist 
in the nucleus of luminous infrared galaxies (the so called LIGs which are characterized 
by extreme mid- to far-IR luminosities L > 1O 11 L and high brightness temperature radio 
emission). The question to assess the relative contribution of each form of activity to the 
LIGs energetics, and to determine a potential evolutionary progression from one to the 
other has been recently addressed by Smith, Lonsdale & Lonsdale (1998, hereafter SLL98). 
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Searching for a SB-compact object AGN connection in these luminous infrared galaxies 
they have demonstrated that nuclear starbursts can in principle explain the observed VLBI- 
and VLA- scale structures and the high infrared luminosities of most of these objects, 
although in some cases severe constraints must be imposed on the initial mass function 
and, in general, complexes of extremely luminous radio supernovae (RSNs) (clustered in 
pc-scale clumps) are required to explain the high-Tf, observed. In a sample of 10 LIGs, they 
found potentially 6 SB systems, 3 possible hybrid/intermediate systems and 1 quasar-like 
object (or BH candidate). This would favor an evolutionary interpretation in which a SB 
eventually begins to fuel a coalescing central compact object. Galactic mergers are believed 
to be responsible for the origin of LIGs, and numerical simulations (Mihos & Hernquist 
1996) suggest that the SB occurs relatively late in the merging process (~ 10 8 yr). SLL98 
results, in turn, suggest that a SB phase might carry on for another ~ 10 s yr as the nuclei 
coalesce and a compact AGN develops. Once the compact AGN fully turns on, it must 
quickly consume its shroud of dust thus extinguishing the star formation process. 

A prototype of the LIG class, Arp 220, seems to be dominated by a massive nuclear 
starburst imbedded in a rotating disk-like structure of dense molecular gas with a radius 
~ 235 pc, and with densities (n ~ 10 4 cm~ 3 ) and temperatures similar to those of the 
GMC core regions in our Galaxy (Scoville, Yun, & Bryant 1997, hereafter SYB97). Several 
other systems (e.g., NGC 6240 and Mrk231) bear striking similarities to Arp 220, which are 
suggestive of a centrifugally supported disk of dense molecular gas. These disks may play 
a critical role in the evolution of ultra-luminous galactic nuclei and in their energy release. 
The high density of the gas will undoubtedly promote extremely high rates of massive star 
formation since the conditions are very similar to those of GMC cores. In these cases, 
the presence of strong, highly ordered magnetic fields could provide both the conditions 
appropriate for RSN formation and production of collimated outflows. 
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Nuclear outflows have been observed in several Seyfert galaxies and LIGs (e.g., Arp 
220, NGC 6240, NGC 2782, NGC 4945, and NGC 4235) (see, e.g., Heckman, Lehrnet, 
& Armus 1993, Colbert et al. 1996, Jogee et al. 1998, and references therein). In the 
case of the LIGs, the detection of emission line profiles with typical line widths of 200-600 
km s _1 , and in some cases with FWHM as large as 1000 km s^ 1 ( e.g. NGC 6240), are 
signatures of the presence of high-speed, supersonic outflows where the lines are produced 
by shock-heated ambient gas which is swept by the outflow (Heckman, Armus, & Miley 
1990, herareafter HAM90). In the case of the Seyfert galaxies, radio observations of a 
distance-limited sample of edge-on galaxies by Colbert et al. (1996), have shown that 6 of 
the 10 observed galaxies, have extended radio structures with radial extent between 5 and 
30 kpc. They are oriented at skewed angles (A = 35° — 90°) with respect to the galaxy 
minor axis and do not look like the wide-angled winds observed in typical starburst galaxies. 
Their morphologies resemble instead, those of more collimated outflows originating in the 
nuclear region (within 1 kpc), like the jets observed in extended radio galaxies, although 
the Seyfert outflows are weaker (with radio power at 4.9 GHz, P 49 = 10 20-22 W Hz -1 ), 
smaller, and less collimated than those of the radio galaxies (which have typical radio power 
at 1.4 GHz, P L4 > 10 24 W Hz" 1 ). Terminal outflow velocities of ~ 25 to 150 km s" 1 and 
magnetic fields ~ 10~ 7 to 10~ 6 G at 5 to 30 kpc scales (implying average fields ~ 10~ 5 G at 
~ 1 kpc) have been inferred from the observations. 

Based on the evidences above, it seems plausible that the LIGs and perhaps some of 
the low-level-activity AGN class of Seyfert galaxies may be powered by nuclear SBs at least 
during part of their active lives. In the present work, we assume that a SB, which has been 
induced in an interaction or merging process, is in progress in the nuclear region of these 
objects and examine the conditions under which collimated outflows can be produced in 
such a scenario. Previous purely hydrodynamic numerical studies of starbursts generated 
within few kpc radius from the nucleus have shown that SB mass loss dynamics may 
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provide favorable conditions for a free vertically directed escape of hot gas from the galactic 
environment at the vigorous SN dominated phase of the SB (e.g., Tomisaka & Ikeuchi 1988, 
De Young & Gallagher 1990, Suchkov et al. 1994, Tenorio-Tagle & Muhoz-Tuhon 1997). 
Those models have successfully explained the formation of the poor-collimated, wide-angled 
super- winds (with radius ~ 1 kpc) often observed in SB galaxies (e.g., HAM90). In 
this work, however, we are mainly interested in systems with potential SB-AGN activity 
connection and thus examine the production of outflows within radial distances of few 
tens of parsecs from the central core. We consider a mechanism that can provide efficient 
collimation and acceleration of the emerging outflow on these scales. We assume that part 
of the gas released by star formation in the rotating SB system is quickly collected into a 
disk with similar conditions, for example, to those observed in the Arp 220 nuclear disk. 
Observations support the evidence for magnetized outflows. Since magnetic winds can be 
produced by rotating objects with a magnetic field anchored into them (e.g., Weber and 
Davis 1967), we assume that magnetic fields are present due to magnetic-flux captured 
from the environment during the formation of the disk, and compression by the accretion 
process (see, e.g., Spruit 1996). 

The perception that rapid rotation coupled with strong magnetic fields could drive 
outflows was first realized in the context of the solar wind theory (e.g., Weber and 
Davis 1967, Mestel 1968). Later, Blandford & Payne (1982) introduced a scale free 
hydromagnetic wind model for jet acceleration from a magnetized disk which gives 
centrifugally driven outflow solutions. These solutions also have the possibility of being 
self-collimated (Heyvaerts & Norman 1989). (See, however, Okamoto 1998 for a discussion 
on the limitations of these models.) Similar outflow models were subsequently applied to 
proto-stellar objects (e.g., Pudritz & Norman 1983, Pelletier & Pudritz 1992, Wardle & 
Konigl 1993, Shu et al. 1994, see also Spruit 1996, for a review), galactic X-ray binaries 
(e.g., Blandford 1993), and proto-jovian planets (Quillen & Trillen 1998). In the present 
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work, we use the Pelletier and Pudritz solutions (hereafter PP92) to investigate the 
formation of magneto-centrifugally accelerated outflows in a SB environment. 

In §2 we try to build a consistent SB scenario and examine whether massive nuclear 
star formation may be present at high enough rates to drive galactic flows out from the 
surface of the nuclear gas disk produced in the course of the SB evolution. In §3, we discuss 
the physical conditions for the nuclear disk to produce collimated outflows based on the 
observed characteristics of the Arp 220 nuclear disk, and estimate the radius, accretion 
and wind outflow rates, terminal velocities, kinetic energy rates, and magnetic field ranges 
for which magneto-centrifugal outfows could operate. In §4, we draw our conclusions and 
briefly discuss our results. 

2. A Starburst Scenario 

In an attempt to build consistent starburst scenarios and evaluate their energetics, 
mass loss, and rate of supernova production, in previous work (de Gouveia Dal Pino & 
Medina Tanco 1997, hereafter GT97) we have calculated the evolution of several starburst 
stellar clusters based on the evolutionary population synthesis models of Schaller et al. 
(1992). As initial conditions we have assumed a Salpeter initial stellar mass distribution 
function [$(m) oc m~ 7 ], with a spectral index 7 = 2.35, and stellar masses in the range 
[Mi = 1 - 8 M , M u = 120 M ]; initial number of stars N c = 10 6 -4 x 10 9 (corresponding to 
cluster radii R c ~ 0.2 — 705 pc and masses M c ~ 10 7 - 10 10 M Q ); and an energy deposited 
by each supernova (SN) E = 10 51 erg s -1 . Assuming also that the system was rotating 
with an ellipticity e = J/r c v s ~ 2.5 x 10" 3 , where J is the angular momentum per unit 
mass (Salpeter and Saslaw 1996, Begelman & Rees 1978), we have found that part of the 
gas released by star formation may cool quickly by free-free emission in time scales much 
smaller than the evolutionary time scales of the system (r <^ 10 s yr) and thus fall towards 
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the centre and settle into a disk with an initial characteristic radius r d ~ e 1 / 2 r c . In fact, 
rapid fueling is expected by the large star formation rates and gas conversion into stars 
detected in nuclear SB regions (Kennicutt et al. 1987, Lo et al. 1987). These observations 
indicate that fueling should occur in time scales smaller than the conversion time into stars 
(~ 10 7 — 10 8 yr) if the gas is to accumulate faster than it is consumed (e.g., Tenorio-Tagle 
& Munoz-Tuhon 1997). As an example, Figure 1 (see also GT97) shows the results for a 
cluster with 3 x 10 8 M and a core disk which eventually becomes also gravitationally 
unstable and undergoes star formation (see below). For a cluster-disk system with masses 
and radii about ten times larger (smaller) than those of Fig. 1, all the luminosity values 
and supernova rates increase (decrease) by an order of magnitude. 

As investigated in previous SB-driven outflow models, a strong nuclear SB can, in 
principle, drive a wide-angled wind out of the system (see e.g., HAM90). We will first 
discuss this process in terms of the energy available to produce winds in a SB system and 
then, on subsequent sections, examine the mechanism for providing outflow acceleration 
and collimation. If an outflow is driven by a SB, then the observed non-thermal radio power 
{P n t) should be correlated with the supernovae production rate (i>sn) in the SB through 
(Condon & Yin 1990): 



where v is the observed radio frequency, and a r ~ 0.7 — 0.8 is the spectral index. Using 
the observed radio power for the Seyfert outflows we find (see Condon et al. 1996) 
v SN ~ 5.1 x (10~ 3 — 10 _1 ) yr -1 , which is consistent with the predicted values of u SN by 
the SB models (see, e.g., Fig. 1). 

Alternatively, simple analytical prescription can be used to estimate the expected 
supernova production rate vsn as a function of the star formation rate (SFR). For a power 
law initial stellar mass function, with upper and lower stellar mass cutoffs m u ~ 12OM , 




(1) 
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and mi ~ 1M , respectively, we have (e.g., Scoville & Soifer 1991) 




(2) 



where tusn = 8M is the minimum mass for supernova detonation, or 



u SN - 0.5 yr 



-i 



( 



20 M yr 



SFR 



-i 



) 



(3) 



which shows that SFR ~ (5 — 20) M yr 1 result SN rates consistent with those required to 
explain the observed radio power in Sey outflows. 

In the case of the LIGs, the observed infrared luminosities imply star formation rates 
SFR ~ 16 - 125 M Q yr -1 . SBs models with these SFR result in supernova rates u SN ~ 0.1 
- 2 yr -1 (SLL98). These are also consistent with the u$n values required by the observed 
non-thermal radio luminosity in these objects. In the particular case of the Arp 220 source, 
the supernova rate required by the observed non-thermal radio power is u SN ~ 2 yr -1 
(which is in agreement with the predicted rate by a SB model with cluster and disk masses 
an order of magnitude larger than those shown in Fig. 1, Md ~ 1O 9 M ). 

SBs with v S N — 5 x 10~ 3 - 2 yr -1 may inject kinetic energy into the ambient medium 
at a rate (see, e.g., Eq. 1; HAM90) 



which is comparable to both the rate of gravitational energy released by the infalling gas 
and the SN luminosity in the core disk of the SB models (see, e.g., Fig. 1, for which u SN 
^ 0.2 yr -1 , and L S n — L K ~ 10 43 erg s" 1 ), and smaller than (^ 20 %) the observed 
bolometric luminosities of the sources (L bo i ~ L FIR ^ 10 n L o for typical LIGs and Sey 
galaxies). In fact, we will see later (§3) that only a fraction of this SN energy rate is actually 
necessary to power the outflow. The remaining SN energy will be possibly lost in radiation, 
as indicaterd by recent numerical studies (Thornton, Gaudletz, Janka, & Steinmetz 1998) 



L K ~ 3.5 x 10 43 u SN *t (2 x 10 



,41 



7 x 10 43 ) ergs 



i 



(4) 
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3. Modelling the Outflow 



3.1. 



Disk Parameters 



In the previous section, we have assumed that a rotating nuclear disk of gas is produced 
during the SB evolution and have evaluated the amount of kinetic energy that such a 
system is able to inject into the ambient medium. Now, before deriving the conditions for 
outflow acceleration and collimation from the surface of the disk, we must explore fiducial 
physical parameters for the disk. In earlier studies of SB-driven, wide-angled super-winds 
(e.g., HAM90), global physical conditions were derived for the central regions (radius 
1 kpc) of SB galaxies. Here, since we are focusing on the inner regions (located within 
few tens of pes radii from the nucleus) of systems with potential coexistence of a SB-AGN 
activity, we take the observed data of the nuclear disk of the LIG prototype, Arp 220 
(SYB97). High-resolution CO imaging and dust continuum emission of the nuclear region 
of Arp 220 has revealed the presence of a rotating disk feature of molecular gas with mean 
density n H2 — 2 x 10 4 cm -3 , a radius ~ 250 pc, and two stellar nuclei at ~ 235 pc from 
the center. Most of the disk mass resides in the molecular gas, rather than in the stars, and 
most of the molecular gas (IVLj ~ 5.4 x 10 9 M Q ) is concentrated in a very thin disk within 
the double nuclei distance from the center. The fact that the bulk of the molecular gas has 
relaxed into a disk with large mass of gas concentrated interior to the double stellar nuclei is 
consistent with a scenario in which the gas in merging systems settles into the center faster 
than the stellar components of the SB. The mass in the Arp 220 disk is approximately 
distributed according to the empirical relation (see Fig. 8 of SYB97) 



with ~ 1.66 (for r £ 200 pc). 

The rotation velocity is related to the mass distribution through the simple spherical 
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approximation v rot = {GM r jr) 1 l 2 . 

For a self-gravitating disk, the effective thickness (H) of the gas can be obtained from 



90 km s 1 is the vertical dispersion velocity (SYB97), and £ = M d /A d = 2Hp d is the mass 
surface density of the disk. Within a radius r ~ 200 pc and M d ~ 5 x 10 9 M , E ~ 3 x 10 4 
M pc~ 2 for Arp 220, so that 



This implies a thin disk at r ~ 200 pc with H/r ~ 0.075. In the inner region it can be 
thicker however. 

The stability of the gas disk to axi-symmetric perturbations has been examined by 
SYB97 who have used the Toomre (1964) criterion Q ~ ak/^GT^) (where k ~ 2VL is the 
epicyclic frequency). The fact that Q ^ 1 within r ~ 400 pc implies that the gas disk is 
unstable and will develop large density enhancements, ultimately leading to the formation 
of massive star clusters in the disk. This result is consistent with our previous assumption 
of a disk originated in a SB environment also processing star formation. 

The extreme gas densities in the central disk of Arp 220 (n d = n H , 2 ~2x 10 4 cm" 3 ) will 
result in high accretion rates. Assuming a viscous transport of angular momentum, SYB97 
have used the Pringle (1981) treatment to estimate the rate of gas accretion as a function 
of the disk radius. For the typical Arp 220 parameters, they found accretion rates M d ~ 
100 M yr _1 peaking at ~ 250 pc (which corresponds to a kinematic viscosity v ~ 3 x 10 27 
cm 2 s _1 ). 

We now must check whether the above accretion rate M d is consistent with the range 
of temperatures and pressures estimated for the disk. Accretion through the disk implies 
that there is dissipation of energy into the disk. For a self-luminous disk the temperature 



the standard condition of hydrostatic equilibrium (Spitzer 1942), H = ( ^frff ), where 
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as a function of the radius is T r 4 = f|^^§r, where <j S b is the Stefan-Boltzmann constant. 
Using the accretion rate mentioned above we estimate 

Td~ 14 K ( M ' ) V4 ( *■ ) ^ (^-) " m (7) 
d ~ ^5xlO 9 M ; ^lOOM yr-y ^200pc/ W 

which is consistent with the temperatures estimated from the OH observations in Arp 220 
(~ 17-21 K; SYB97). 

The sound speed at this radius would be v s ~ 0.4 (T/14 K) 1 / 2 km s" 1 , and the disk 
pressure 

p ' s4xr "*»1iJS5F)(si)' (8) 

We note, however, that in our SB / disk model, massive star formation may still be occurring 
after the formation of the disk and even in the disk itself. In such a case, one would expect 
the presence of an ionized, warmer gas component produced in interactions between the 
ejected matter in the SB (e.g., HAM90). A disk with a dominant warm gas component 
(with Td ~ 10 4 K) would have a higher pressure Pd — 10~ 8 dy cm~ 2 . In any case, no 
matter whether the inner disk is cold (i.e., presently not undergoing massive star formation) 
or warm, the magneto-centrifugal mechanism for outflow production that will be discussed 
below is applicable to both. 



3.2. Outflow Conditions 

Previous SB-driven modeling of wide- angled super-winds (e.g., Tomisaka & Ikeuchi 
1988, De Young & Gallagher 1990, Suchkov et al. 1994) have shown that the multiple 
action of stellar winds and SN explosions may lead to the formation of a shell of shocked 
disk material which, after blowout, is dragged to the halo, where it forms the walls of a 
bi-conical cavity. This would, ultimately, create conditions for a free directed escape of hot 
gas at the later, more vigorous SN dominated phase of the SB. Formulation derived by Mac 
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Low & McCray (1988) and Mac Low, McCray, & Norman (1989) gives a simple condition 
for a SB-driven outflow to blow out of the atmosphere above the disk. The rate of energy 
released by the SN, Lk (eq. 4), must satisfy 

Lk > 3 X 1O5L0 ( 2 xlOW3) (l^) (lO-*dycm-*) (9) 

which is a condition more than satisfied by the candidate systems, which have 
L K ~ (5.3 x 10 7 - 2 x 10 10 )L© (eq. 4). We note that this condition still holds for more 
compact, denser disks with, for example, H ~ 1 pc, n d £ 10 8 cm -3 , and < 1CT 3 dy 
cm~ 2 . 

The results above suggest that nuclear SBs may satisfy the conditions for driving 
winds out of the galactic core of the Seyferts and LIGs, but, how to explain the collimated 
morphology of the observed outflows that we are investigating and their acceleration? One 
could invoke a thick gas torus surrounding the core to collimate the "wind" in a similar way 
to that proposed in standard models for AGN (see, e.g., Antonucci 1993). As mentioned 
earlier, hydrodynamic models for SB super-winds predict the formation of a cavity that 
provides a natural channel for vertical escape of the gas ejected by the SN. Recently, 
Tenorio-Tagle & Munoz-Tuhon (1997), have explored hydrodynamic solutions in which the 
rapid fall of disk matter towards the center develops a shock front in the inner region of the 
disk, and the material crossing it experiences a splash and is ejected vertically to the lower 
pressure zones above the disk, thus forming a (poorly collimated) kpc-long cone. All these 
models seem to be able to provide appropriate conditions for the creation of the observed, 
poor-collimated super-winds emerging from the kpc-scale galactic disks of SB galaxies. 
Here, in our investigation of the dynamics of nuclear SB/disks, we examine another 
potential mechanism in which efficient collimation and outflow acceleration is provided by 
magnetic fields anchored into the core disk of the SB — an assumption that is supported by 
the observation of magnetic fields in the outflows of the Seyferts (B ~ 3.5 x 1CT 7 — 1.6 x 10~ 6 
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G at few kpc above the disk; Colbert et al. 1996). Such mechanism could also, in principle, 
be operating in the inner (possibly unobservable) regions of wide-angled super-winds of 
normal SB galaxies. 

In a magnetized gas under low ionization fraction condition, for an outflow to be 
launched the neutral component must be well coupled with the ions (e.g. Wardle & Konigl 
1993). This coupling is described by 77, the ratio of the dynamical time to the mean collision 
time of a neutral atom in a diluted sea of ions (e.g. Wardle & Konigl, Shu et al. 1994, 1995), 
77 = where p,i is the density of ions and v rot is the rotational velocity of a particle in 
a circular orbit. When r\ 1 the neutrals move with the ions and so can be described as 
well coupled to the magnetic field. Here 7, ~ 3 x 10 13 cm 3 g _1 s^ 1 is the drag coefficient 
between ions and neutrals and pi, the density of ions, is the product of the atomic mass of 
the ions, the ionization fraction and the total density. We can then estimate 

For clouds with temperatures, abundances and column densities not much different from 
those in the disk examined here, Friaga & Jafelice (1998) have found that the ionization 
state is roughly inversely correlated with the density. From their curves we have estimated 
the corresponding value of n e /riH — 10~ 5 — 1CT 2 for our density and temperature ranges. 
We conclude that over the whole disk r\ ^> 1 and the neutrals can be said to be well coupled 
to the ions. However this does not mean that the particles are necessarily well coupled to 
the magnetic field. 

The magneto-centrifugal wind models can operate where the charged particles are 
sufficiently coupled to the field lines that they are flung out centrifugally from the outer 
layers of the disk (for a favorable field geometry; see below). This occurs when the 
magnetic Reynold's number is Rtu 1. We can estimate the magnetic Reynolds 
number as Rem = Vr °* H , where v M is the magnetic resistivity, v M = c 2 /(A7ca e ), and a e is 
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the electric conductivity determined from collisions between electrons and H 2 molecules, 
°" e — n„ m<o- v> — ll* 10 16 (™e/Vff) s -1 (where a c is the cross section for e — H 2 
collisions, v is the relative velocity between them, and < a c v > is the Maxwellian mean of 
the product; Hayashi 1981, Takata & Stevenson 1996). Using the Arp disk parameters we 
estimate 

*» > » * 1018 {t&sf wy (^) fcr (n) 

so that we can assume that the magnetic field is well coupled to the gas. 

3.3. Magneto-centrifugally Driven Outflows 

Collimated winds may be produced by a rotating disk with a magnetic field (B) 
anchored into it (see, e.g., Spruit 1996 and refereces therein). The magnetic field and disk 
structure will determine how much mass is launched out of the disk. 

As remarked before, we here assume that the field is produced due to magnetic flux 
capture from the environment during the formation of the disk in the SB, followed by 
advection and compression by the accretion process. Poloidal flux capture in this way 
cannot be destroyed by local processes in the disk, it can only escape by diffusing radially 
outward (e.g., Spruit 1996). It is out of the scope of this work to formulate a theory for 
(turbulent) diffusion in the accretion disk, so that we cannot reliably predict the distribution 
of the magnetic field. However, since all the energy densities are expected to increase 
towards the center, it is reasonable to assume that the field will increase inward as well and 
form an opening field geometry. As argued by Spruit, Stehle, & Papaloizou (1995), fields 
with such a configuration can be quite strong, with magnetic energy densities exceeding the 
gas pressure, which must make them ideal for magnetic outflow production. 

Just above the disk, the gas density is typically low enough, so that the magnetic 
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energy density is large compared with the thermal and rotational energies. The field in this 
region may be thus assumed to be force- free (JxB = 0). The poloidal magnetic field lines 
that raise from the disk are somewhat like rigid wires that control the flow. The gas of the 
disk surface will be forced to corrotate with the field lines raising out of the disk. For a 
sufficiently cold disk, the solely forces acting on a parcel of gas will be the gravitational force 
and the centrifugal force due to rotation. Along the field line above the disk, the centrifugal 
force will increase with the distance from the axis and the effective gravity decrease. When 
the centrifugal force component along the line exceeds that of gravity, the gas tied to the 
field line is accelerated outward. This outward centrifugal acceleration continues up to 
the Alfven point [the location where the outflow poloidal speed reaches the Alfven speed, 
v — va — B 2 /(A tv p) 1 / 2 , and stops when the flow attains the fast magnetosonic point. At 
that point, the field is no longer strong enough to enforce corrotation. Beyond the Alfven 
surface (pv 2 > B 2 /8n), the inertia of the gas causes it to lag behind the rotation of the 
field line, and the field winds up thus developing a strong azimuthal component (B^). The 
magnetic stress force associated to the azimuthal field [— B 2 /(4irr), where r is the radial 
distance from the disk axis] will provide the flow collimation (see e.g., PP92, Spruit 1996). 

A self-gravitating accretion disk with an accretion rate Md releases a known amount of 
gravitational binding energy of gas as it spirals inward to the center. The presence of viscous 
stress causes the gravitational binding energy to transform into heating and the angular 
momentum is radially transported outward by turbulent stress. In the presence of magnetic 
fields, however, part of the gravitational energy released in the accretion is converted into 
mechanical energy of a wind. In addition, the angular momentum is efficiently transported 
by the wind through the disk surface, and the wind mass-loss rate M w has only to be a tiny 
fraction of M d (see, e.g., PP92 and references therein). Using the disk angular momentum 
equation one can show that as long as the magnetic pressure is comparable or larger than 
the gas thermal pressure, the magnetic torque will completely dominate the viscous torque 
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(e.g. PP92). In fact, the ratio of the magnetic to the viscous torque is 

i~) = (i^) (^h) (12) 

where ass < 1 is the Shakura-Sunyaev viscous prescription (e.g., Frank, King, & Raine 
1992), and ta is the Alfven radius. Thus, since ta/qssH > 1, and the magnetic pressure in 
the disk (B 2 /8ir ^ 2.5 x 1CT 8 dy cm -2 ; see below) is effectivelly of the order or larger than 
the thermal pressure (Eq. 8), we expect that the angular momentum loss through the wind 
will dominate. 

In magneto-centrifugal disk outflows, the angular momentum of the wind is tied to that 
of the disk, and so the wind flux is directly related to the accretion rate, M d) through the 
disk. To investigate specific conditions under which collimated winds could be generated in 
the SB/disk scenario drawn in the previous paragraphs, we here follow the non-self-similar, 
steady-state, magneto-centrifugally driven wind models developed by PP92. In their 
solutions, the angular velocity of the disk is assumed to be a power law of the magnetic 
flux, Vt oc tp~ a . They have found that for a power law index a = 3, the outflow carries a 
constant current intensity everywhere and is free of the divergence found in self-similar 
models (e.g., Blandford & Payne 1982). We here adopt this relation. Now, according to Eq. 
5, the angular velocity scales with the disk radius as 

Q = { ™f 2 oc r^" 3 )/ 2 oc r" - 67 . (13) 

Therefore, the magnetic flux will also be a power law of the radius, 

iP oc r^ 1 oc r - 223 . (14) 

This, in turn, implies the following relation for the poloidal magnetic field in the disk as a 
function of the radius (PP92), B r oc r~ s , with 5 = 2, or 

B r = B M (^y (15) 
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where ~ 10 5 G is the average value inferred from the observations of the Sey outflows 
at the kpc-scale, = 1000 pc. 

In these models, the ratio of the radius r to the Alfven radius, ta, ^(r/r A ) 3 , can be 
determined by the parameter (PP92) 



M d VGW r 



r 



For a wind to be produced the field lines must be more than 20° from normal to the disk 
and e < 1. Assuming M d constant, the relation above for B r (Eq. 15) and Eq. 5 for the 
disk mass result 

e oc r^+ 25 oc r 2 - 33 . (17) 

Using these relations and the physical conditions in the Arp 220 disk (with constant 
M d = 100 M Q yr' 1 ) , we can determine a maximum radius r e for which e(r e ) = 1, 

M \ -0 ' 43 / R \ - 86 / r \ 1 - 72 



This value determines a characteristic maximum radius for the outflow zone in the disk. 
We note, however, that it is sensitive to the adopted value of the magnetic field on the large 
scales Boo which is only barely known. (Notice that Eq. 15 implies B(r e ) ~6x 10~ 4 G on 
the disk surface at r e .) 

For these models the mass flux in the wind, M W) is related to the accretion rate through 
the disk by (PP92) 

M w = fM d (19) 

where 

2 



(20) 



HAM90 have estimated values for M w for wide-angled super-winds observed in 
far- infrared galaxies ranging from a few M yr _1 in the low-power sources (e.g., M82), up 
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to ~ 50 M yr _1 in powerful sources (like Arp 220 and Mrk 273). These values correspond 
to the total amount of mass loss rate integrated over the whole radius of those winds. In 
the inner, few tens of pc scales, M w is expected to be smaller. In fact, in our model, M w 
approximately scales as M w oc ^ 2 (eq. 6.7 of PP92), or according to Eq. 14 above, M w 
oc r 045 . Therefore, the value M w ~ 50 M Q yr -1 for the Arp 220 outflow will correspond to 
disk radius >> 100 pc . We can then assume that in the inner few pes M w ~ 10 M yr _1 
and take M w /M d ~ 0.1. This value is compatible with estimates for YSOs and AGN outflow 
models (see, e.g. PP92 and references therein). Substituting this value into Eqs. 19, and 20 
above, we can then estimate a fiducial value for the ratio r A /ri ~ (M w / M d )~~ 1 / 2 ~ 3.2 at a 
fiducial radius r\, and according to Eq. 16, 

€i ~ ^(^) 3 - 1-1 x 10- 2 (21) 
Thus substituting this value into Eq. 17, 

/e n 0.43 

r i = r e(f) - 18pc, (22) 

or (according to Eq. 16) 

c ~ 1 1 x ID- ( Md ) ( Mr > ) V2 ( Br * ) " ^ (23) 

VlOOM yr-V V9.2 x 1O 7 M / V3.1xlO- 2 G/ V18pc/ ' 1 1 

The terminal speed of the outflow can be evaluated from the solution of the Bernoulli 
equation for the flow in the asymptotic domain (e.g., PP92, eq. 6.6) 

Voo {r) ^ (fe) r 2 oc r" - 45 (24) 

<J e { r e) 

where (f2 e r e ) is the rotation velocity in the disk at the outer outflow radius r e , and A is a 
constant everywhere in the wind that can be estimated from Eq. 6.7 of PP92, and Eq. 14 
above 

~ 9.35 (25) 



A 



4^ 2 M d 
2 e[r ' e) MJr e ) - MjrA 



0.45 
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where we have used the approximate scaling law for M w (r) oc ip 2 oc r a45 . Substituting A 
into Eq. 24, we find the expected terminal velocity at the radius r\ at which most of the 
wind originates 

^- 973kmS " 1 (28&)(^j)"°' 45 < 26 ' 
and t>oo(r e ) ~ 0.42 foo(rj) ~ 408 km s -1 . These estimates are consistent with the detection 
of emission lines typically with widths ~ 200 — 600 km s -1 in LIGs. They indicate the 
presence of a high speed pre-shocked wind gas that was possibly decelerated to these 
velocities. In the case of Arp 220, the observation of line widths declining from ~ 700 km 
s _1 near the axis of the outflow to ~ 300 km s^ 1 at the outer radius are also consistent with 
the prediction of our model of a decreasing outflow velocity with radius. On the other hand, 
the estimates above are higher than the observed outflow speeds in the Seyfert outflows 
(voo, obs— 25 — 150 km s _1 ; Condon et al. 1996). Nonetheless, according to Eq. 24, these 
observations can be fitted by a more compact disk (like that of Fig. 1, for instance) with 
a somewhat larger ratio M w /Md (Eq. 21). Eqs. 21, 22, 24 and 25 indicate that disks with 
r e <^ 10 pc and M w {r-i) / Md ^0.1 will provide Voo in better agreement with the observed 
Sey outflows. 

In summary, the results above indicate that the observed terminal outflow velocities 
are consistent with an outflow generation in the described SB-disk scenario with wind mass 
loss to disk accretion rate ratios ^0.1. 

One of the attractions of a magnetically driven outflow model is that it is able not only 
to produce and collimate the outflow, but also to efficiently extract angular momentum 
from the disk, allowing it to collapse (e.g. Blandford 1976, PP92). Just a small fraction of 
the infalling mass, like that estimated above, is required to be carried out by the outflow 
for it to take out most of the angular momentum of the system. Further, in the outflow 
model investigated here, it is easy to show that the total kinetic power carried by the 
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wind is exactly one-half of the gravitational energy released in the accretion disk, i.e., 

t _ lr _ 1 GM re M d T- /- n 0.451 
■L'w 2 acc 2 r e 



1 _ li 



, where M re is the disk mass within r e . Again, using 
the Arp 220 parametrization, we find that 



/ / \ 0-45' 

1 _ ( n/r e \ 



0.14 y 

(27) 

This can be compared with the SN energy rate estimated for the Arp 220 disk, L K ~ 7 
x 10 43 erg s _1 (Eq. 4), which results L W /L K ~ 2 x 10 -2 , so that the total kinetic power 
carried by the magneto-centrifugally driven outflow is much smaller than the available SN 
energy rate predicted by the SB models. 

The fact that the open field lines must have an inclination angle with respect to the 
disk not larger than 9 = 70° (PP92) results that the height za of the Alfven surface is 
z a < r A tan 9, or 

* < ^fr*) fe) ■ (28) 

We can also estimate the outflow number density at the Alfven surface (and beyond) 

3-/3 

from the condition that n A oc ip 2 (PP92), and according to Eq. 14, ua oc r~ oc r 0A5 . 
Using the number density at the base of the disk, which was estimated from the Arp 220 
disk, n{r-i) ~ 3 x 10 4 cm~ 3 , the relation above implies that 

/ I \ °-45 / \ 0.45 

„^ 5xl 0^ m -(^) (^) . (29) 

According to PP92 models, the flow will recollimate beyond the fast-magnetosonic 
(fm) point. The Mach number at this point is 

M fm = (2 3 / 2 A 1 / 2 - 3) 1 / 2 (30) 
Substituting the value of A given by Eq. (25) we find M/ m ~ 2.4. Thus flows achieving 
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modest fm-Mach number ~ 2 must recollimate and form magnetic focal points beyond 
about a hight za ~ 150 pc. 

4. Conclusions and Discussion 

SBs with a disk core formed by infalling gas produced during star formation in a 
rotating stellar cluster, seem to be able to provide the amount of energy required to drive the 
weak, directed outflows of Seyferts galaxies and LIGs. In the latter source class, for example, 
there seems to be strong observational support for the presence of nuclear SBs (e.g., SLL98). 
We have found that massive nuclear SBs with cluster masses M c ~ 10 9 — 1O 1O M , disk 
masses M d ~ 10 8 — 1O 9 M , and supernova rates u S n — 5 x 1CT 3 — 2 yr" 1 (which are 
consistent with the inferred i>sn from the observed non-thermal radio power) may inject 
kinetic energies ~ 10 41 — 10 43 erg s -1 which are high enough to blow out directed winds 
from the accreting disk surface (§2), within the SB lifetimes. 

In our models, the acceleration and collimation of the outflow is provided by moderate 
magnetic fields anchored into the rotating SB-disk. Based on empirical conditions derived 
from the disk of the prototype LIG, Arp 220, and the magnetized outflows of Sey galaxies 
and LIGs, we have found that magnetic fields with intensities ^8x 10~ 4 — 5 x 10~ 2 G 
are able to accelerate the flow up to the observed terminal velocities (^ few 100 km s -1 in 
the case of the Seyferts, and ~ 400 — 950 km s _1 in the case of the LIGs), with wind mass 
loss to disk accretion rate ratios M w /Md ^ 0.1 (where M d ~ 100 M Q yr -1 ). The outflow 
carries only a small fraction (few percent) of the available SN kinetic energy rate, and is 
produced within a wind zone in the disk of radius ^ 100 pc in the case of the LIGs, and 
^ 10 pc in the case of the Seyfert outflows (§3). The model also predicts a decay of the 
terminal velocity with radius which is compatible with observations. 
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The formation of an outflow in the proposed SB-magnetized disk scenario depends on 
a series of conditions which are evolutionary linked, such as, the initial amount of angular 
momentum and ellipticity in the central stellar cluster that will provide the formation of 
the rotating disk of gas from the SB; the development of an appropriate magnetic field 
geometry with open lines during disk collaps; and the ionization state of the gas which must 
be large enough to couple the magnetic field to the gas. SB-driven outflows are observed 
not only in IR luminous nuclear SB systems, but also in other galaxy types, as for example, 
blue compacts and irregular dwarf galaxies like NGC 1569 (e.g., Greggio et al. 1998). 
The collimation-acceleration mechanism here described could, in principle, also apply to 
the inner regions of the wide-angled super-winds observed in these systems. They could 
also anchor a nuclear magnetized disk generated during the SB if the initial conditions 
for the SB were not very much different from those described here. However, in the case 
of the amorphous dwarf systems like NGC 1569, which is actually a post-SB system, it 
is not unlikeley that the SB-driven large-scale wind had stripped the dwarf of most of 
its gas before the gas had time to partially collect into the center and develop a nuclear 
rotating disk and a more collimated, coherent outflow portion. In forthcoming work, we 
aim to explore self-consistent numerical models to examine the dynamical evolution of 
outflows produced by a rotating magnetized disk originated in different SB environments 
and determine their range of applicability. For that purpose, we will employ a modified 
version of a 3-D numerical SPH code originally designed for supersonic jet studies (e.g., de 
Gouveia Dal Pino & Benz 1993, de Gouveia Dal Pino & Birkinshaw 1996, Cerqueira & de 
Gouveia Dal Pino 1998). 

Pure hydrodynamic models (e.g., Suchkov et al. 1994) make specific predictions for the 
observed abundance pattern in the SB-driven super-winds. In those models, the material 
at the wall of the vertical cavity, for example, is mostly composed of unprocessed ISM gas 
that has been shock-heated by the emerging outflow, while the enriched material that is 
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ejected by the SN should fill the interior of the cavity. Similarly, as a collimated magnetized 
outflow impinges supersonically the ambient medium, it must create a bow shock pattern 
at the outer edge (see, e.g., Cerqueira, de Gouveia Dal Pino, & Herant 1997) where most 
of the shocked material that produces the emission lines is made of ISM gas with no 
significant abundance enrichment. The observed optical emission lines should then trace 
normal abundances just like predicted in hydrodynamical models. The hotter, enriched 
material that is ejected by the SN should emerge mainly along with the collimated outflow. 

We note that in our disk model the magnetic pressure in the disk may exceed the gas 
pressure. This may cause the disk to become unstable to magnetic buoyancy. In such a case, 
the magnetic field lines may arise in loops from the disk surface and undergo reconnection. 
This could in turn provide the formation of hot filamentary and clumpy structure on the 
disk corona. It is interesting to note that a patchy disk corona is required by the X-ray and 
7-ray observations of the nuclear region of Sey galaxies in order to explain their complex 
high-energy spectra (e.g., Johnson et al. 1997). 

Finally, we should inquire about the fate of the SB disk (see also GT97). We have 
seen that it may rapidly become self-gravitating (§2.2). As a consequence, it is expected 
to suffer fragmentation and star formation. Either, the disk may become bar unstable, 
a condition which is met if 5b = T rot /W > 0.14, where T rot is the rotational energy and 
W is the gravitational potential energy of the disk (e.g. Shlosman et al. 1989). Our SB 
disks result 5b = 0.5, so that they can, in principle, become unstable to bar formation and 
undergo fission into two components which eventually may collapse to a supermassive BH 
(Begelman, Blandford, & Rees 1980). Based on the fact that the formation of BHs seems to 
be plausible in a SB, we can speculate that both, the standard BH and the SB mechanisms 
may be occurring in different stages of the AGN evolution. In this evolutionary scenario, 
the stronger more collimated jets of radio loud QSOs, for example, could be produced in 
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the higher-level activity phase of the BH accreting the remains of the disk material. This 
hypothesis is examined elsewhere. 
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Figure Caption 

Figure 1. a) Evolution of a stellar cluster and its core disk, with a Salpeter 
IMF (Mi = 1M Q , M u = 12OM ), iV c = 3 x 10 8 , r c = 53 pc, total cluster mass 
M c = 9.5 x 1O 8 M , R d = 2.6 pc, total number of stars in the disk N d = 9 x 10 7 , and total 
stellar mass in the disk M d = 1.6 x 10 8 M Q . Luminosity sources shown: dotted - supernovae 
in the cluster and disk; solid - rate of energy released by the infalling gas in the disk; dashed 
- starlight in the cluster and disk; and dot-dashed - the Eddington luminosity limit in the 
disk, b) The corresponding rates of supernovae production vsn i n the cluster and the disk. 
For a system with masses and radii about ten times larger (smaller) than those above, the 
luminosity values and supernova rates increase (decrease) by an order of magnitude. 



